New Measurement of the Cosmic-Ray Positron Fraction from 5 to 15 GeV 
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We present a new measurement of the cosmic-ray positron fraction at energies between 5 and 
15 GeV with the balloon-borne HEAT-pbar instrument in the spring of 2000. The data presented 
here are compatible with our previous measurements, obtained with a different instrument. The 
combined data from the three HEAT flights indicate a small positron flux of non-standard origin 
above 5 GeV. We compare the new measurement with earlier data obtained with the HEAT-e* 
instrument, during the opposite epoch of the solar cycle, and conclude that our measurements do 
not support predictions of charge sign dependent solar modulation of the positron abundance at 5 
GeV. 

PACS numbers: ll.30.Pb, 95.35.+d,96.40.Kk,98.70.Sa 



Space-borne and high altitude balloon experiments 
have collected a considerable amount of data on the all- 
electron component in the cosmic radiation. Between 
several MeV and about 50 GeV the energy spectra of 
electrons and positrons have been observed separately 
with sufficient precision to permit comparison with pro- 
duction and propagation models. 

The relative abundances of positive and negative elec- 
trons indicate that the majority of the electron com- 
ponent (e*) consists of negative electrons. They are 
thought to be accelerated in the same Galactic sources 
that also generate the nuclear cosmic rays, but their 
observed energy spectrum is strongly affected by syn- 
chrotron and inverse Compton energy losses during prop- 
agation. 

An additional component of electrons and positrons 
in nearly equal proportion amounts to about 10% of 
the total electron intensity, and is attributed to the de- 
cay of secondary particles (mostly pions) generated in 
hadronic interactions of cosmic-ray nuclei in the inter- 
stellar medium. These positrons constitute only a small 
fraction (< 0.5%) of the total observed cosmic-ray inten- 
sity, yet if they are purely secondaries they can be used 
as an effective probe of cosmic-ray propagation through 
the Galaxy. Their fraction is then expected to decline 
slowly with energy because of the declining path length 
of the primary nuclei at high rigidities. 

Recent observations 1, iSISl confirm the almost ex- 
clusively secondary nature of positrons up to a few GeV. 
However, a possible structure in the positron fraction 
near 8 GeV has been observed with the HEAT-e^ in- 
strument 0, 1^ which may defy a simple explanation. 

It has been suggested that a small positron component 
could originate frornparticle interactions in nearby astro- 
physical sources |a, LD, la ISl or may be generated through 



the annihilation of dark matter particles in the Galac- 
tic halo HSIilElElElElEiElElEl such 

a primary positron component could lead to observable 
features such as those indicated by the HEAT-e^ mea- 
surements. 

A new version of the HEAT instrument, HEAT-pbar 
was designed to observe the high-energy cosmic-ray an- 
tiproton flux but it is also suited for the observation of 
electrons and positrons at energies below « 15 GeV. The 
instrument utilizes a multiple dE/dx vs. rigidity tech- 
nique to identify cosmic-ray particles by mass and charge. 
The ionization loss of relativistic particles is sampled in 
two stacks, each of 70 multiwire proportional chambers, 
filled with a Xe/CH4 mixture. They are located above 
and below a central superconducting magnet spectrom- 
eter which measures the particle's rigidity and charge 
sign ji^]. Scintillators at the top and bottom of the de- 
tector system measure the time of flight and, together 
with a scintillator just above the spectrometer, form the 
trigger. Unambiguous discrimination between protons, 
antiprotons, positrons, electrons, tt"*" and /.t"*", tt" and 
/i^ , is achieved |2l| . This instrument was launched from 
Ft. Sumner, NM on June 3, 2000 and was at float al- 
titude for 22 h, at an average atmospheric overburden 
of 7.2 g/cm^. The vertical geomagnetic rigidity cut-off 
along the flight path varied little and averaged 4.2 GV. 

The time of flight (ToF) system measures the particle 
velocity P — v/c with a resolution of — 0.09, permit- 
ting complete rejection of upward-going particles. The 
scintillator signals measure the magnitude of the parti- 
cle's electric charge and select singly charged particles, 
with a resolution of az = 0.14 (in charge units) for each 
counter. 

The particle's sign of charge and rigidity R=pc/Ze are 
determined with the magnet spectrometer which has a 
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field of about one Tesla. We retain only events with at 
least 13 usable tracking points in the magnet's bending 
plane (out of a total of 17) and at least 6 points in the 
non-bending plane (out of a total of 8). 

Only the smaller 50% of the ionization signals in the 
140 chambers of the multiple dE/dx system are retained 
to form the restricted average signal {dE / dxres) ■ Ow- 
ing to higher electromagnetic energy losses, e^ induced 
events tend to result in a greater number of hits in the 
ionization chambers of the dE/dx system than heavier 
singly charged particles. Therefore we only select events 
with a large total number (64 or more) of channels used in 
calculating {dE / dxres) and we require a minimum num- 
ber of hits in the entire dE/dx stack (155). 



120 



100 - 



> 



3 

S5 



30 40 50 60 70 80 90 100 110 
<dE/dx^gj.> [arbitrary units] 

FIG. 1: Distribution of the average restricted energy loss for 
particles in the rigidity range 4.5 - 6.0 GV identified as neg- 
atively charged. The Gaussian functions were fitted prior to 
the final event selection in order to obtain a well-defined as- 
sociation between average restricted energy loss and particle 
species. After the final event selection is applied, non-electron 
events are significantly suppressed (lower histogram). AH par- 
ticles within the shaded area of the histogram, determined 
as described in the text, are selected as electrons (and as 
positrons in the corresponding positive rigidity bin). 

Distributions of {dE / dxres) for particle tracks having 
satisfied all of the criteria are produced in three rigidity 
bands (4.5 - 6.0, 6.0 - 8.9, and 8.9 - 14.8 GV). To illus- 
trate the particle separation power of the instrument, the 
{dE / dxres) distribution for negative particles with rigidi- 
ties between 4.5 GV and 6.0 GV is shown in Fig.^ Elec- 
trons are clearly separated from muons and pions whose 
small mass difference can not be distinguished by our 
instrument. Each of the three negative {dE / dxres) dis- 
tributions is fitted with a sum of Gaussian functions and 
the mean and standard deviation for the 7r~//i~ mass 
peak (m^^, (T,r/i) and the electron mass peak (me, a^) 
are determined in each of the three rigidity bins. We 
define a particle as an electron (or positron) if its re- 
stricted average dE/dx signal falls between (m7rp-)-3(T7r/j) 



and (me-t-ScTe) in the respective rigidity bin. To further 
improve the discrimination against protons and mesons 
we group the 140 proportional chambers into ten mod- 
ules and require that the average measured energy losses 
in each module exceed a given threshold. This selection 
strongly suppresses protons and mesons while it retains 
almost all electrons (see Fig. QJ. The residual contami- 
nation due to proton/muon/pion spillover was obtained 
by fitting and extrapolating the p//u/7r distributions in 
Fig. [3 The background thus estimated for the three 
rigidity bins are 0.4, 0.1, 1.0 events, respectively. Note 
that the selection criteria employed here were tuned to 
optimize the identification and statistical significance of 
the e^ peaks, and therefore differ from those employed 
in the antiproton analysis of the same data set 
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FIG. 2: Distributions of average restricted energy loss for 
events after all selection criteria have been applied. The three 
rigidity bands, 4.5 - 6.0 GV, 6.0 - 8.9 GV, and 8.9 - 14.8 GV, 
are shown from top to bottom, respectively. We accept events 
under the shaded area as positrons. The strong relativistic 
rise in the energy loss for hadrons and mesons compared to 
electrons (which are already heavily relativistic at GeV ener- 
gies) results in the mass peak of the particle species moving 
closer together with increasing rigidity and ultimately limit- 
ing the particle identification. 

The atmospheric overburden during the HEAT-pbar 
flight varied between 4.5 and 11 g/cm^, for an average 
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overburden of 7.2 g/cm^, which resulted in significant 
numbers of atmospheric secondaries. A Monte Carlo sim- 
ulation, described elsewhere l^^, is used to obtain the 
correction factors for atmospheric positron and electron 
production. The corrections for the positrons vary be- 
tween 44% and 52%, depending on energy, and those 
for the electrons between 5% and 6.3%. Uncertainties 
in these atmospheric corrections result in a systematic 
uncertainty of about ±0.01 on the positron fraction. 

The raw e^ particle counts are obtained by tallying 
the events with a restricted average dE/dx signal that 
falls within a selected region as described above and in- 
dicated in Fig. 121 These counts are summarized in Table 
^(Ne+ and Ng- .) They are corrected for the atmospheric 
secondaries, and the positron fractions calculated. 

Finally, a correction factor of 1.109 is applied to rigidi- 
ties at the instrument, to account on average for radiative 
energy losses by the electrons and positrons, and to cor- 
rect to the top of the atmosphere. The correction is done 
in an average way, rather than on an event-by-event basis, 
because the events are grouped in narrow rigidity bins at 
the instrument prior to their identification as electrons 
or positrons based on their energy losses. This factor is 
calculated as a*/"'"^, where t is the average atmo- 
spheric depth in radiation lengths and a — 3.1 is taken 
for the spectral index of the primary electron flux. 

Fig. |31 shows the final corrected positron fractions, 
e+Z (e+ -f e~), as a function of energy, compared with the 
previous combined HEAT-e* results |3| from the 1994 
and 1995 flights. The accessible energy range for our 
measurements is limited towards the lower energies by 
the geomagnetic cut-off and at high energies by the par- 
ticle separation capability of the instrument. The result 
presented here is in agreement with the HEAT-e^ mea- 
surements obtained with a very different technique. The 
result is also compatible with the predicted ratio from 
[2^ which assumes that all positrons are of secondary 
origin (dashed curve in Fig. O . Although the new data 
set weakens the case for a feature at w 8 GeV, an excess 
of positrons in this energy range cannot be ruled out. 
The combined dataset of the HEAT-e=^ and HEAT-pbar 
experiments is shown in Fig. ^ along with a compilation 
of recent measurements and predicted flux ratios. 



TABLE I: Summary of event selection results and the calcu- 
lated positron fraction (in 10~^). EtoA is the particle kinetic 
energy at the top of the atmosphere. Ng+ and N^- are the 
number of observed positrons and electrons for each energy 
bin, respectively. N^°'' and N^l'' are the extrapolated number 
of positrons and electrons at the top of the atmosphere. 
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FIG. 3: The positron fraction measured by HEAT-pbar com- 
pared to measurements with a different instrument (HEAT-e* 
U); The curve shows the expectation from the model calcu- 
lation in (23l . 

We compare our data to a computation of the cosmic- 
ray secondary positron (and electron) spectrum in a diffu- 
sive model for Galactic cosmic-ray propagation provided 
by [22] for different source injection spectra. With in- 
jection spectra based on the locally measured nuclear 
cosmic-ray spectra, the data follow the general trend of 
the model prediction up to a few GeV, beyond which the 
observed positron fraction is higher than the calculated 
one. This enhancement is plausible in a variety of differ- 
ent theoretical scenarios. 

Motivated by a measured diffuse gamma-ray spectrum 
above 1 GeV that is much harder than expected jl^l, lo- 
cal cosmic-ray fluxes have been calculated based on hard 
proton injection spectra 23]. Interestingly, this does not 
only lead to a better flt of the diffuse gamma-ray data but 
also to an enhancement in the positron fraction above a 
few GeV. However, measurements of cosmic-ray antipro- 
tons provide rather substantial evidence against the idea 
of explaining the diffuse gamma-ray excess with a hard 
nucleon spectrum Therefore, if future measurements 
confirm the positron excess, then there may be additional 
sources for the cosmic-ray positrons. 

An interesting suggestion is that annihilating dark 
matter particles in the Galactic halo, perhaps supersym- 
metric particles, are a source of high-energy positrons [lol 
[nlEllllElESElEElEi.and references therein) . 
In conventional models, the lightest supersymmetric par- 
ticle (e.g., neutralino) is assumed to have a mass larger 
than the W mass and direct decay will produce a sub- 
stantial positron signal at about half the neutralino mass. 
Subsequent heavy lepton and b decays could result in an 
enhancement of the positron fraction at much lower en- 
ergies, as shown in Fig. HI [l5|. 

Other possible contributions to the cosmic-ray positron 
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FIG. 4: The positron fraction as a function of energy 
for the combined HEAT-e* and HEAT-pbar data , com- 
pared to model predictions and other recent measurements 
( CAPRICE 0, Golden AMS H). Dates in parentheses 
give the year of the measurement and not the publication. 
The solid curve is the positron fraction based on a purely sec- 
ondary production of positrons given by . The dashed and 
dot dashed curves are the ratios including contributions from 
Higgsino LSP decay 15| and gamma-ray pulsars m, respec- 
tively. 



flux have been proposed, including synchrotron produced 
e^ pairs from Galactic pulsars 6]. Whether the appar- 
ent excess in the HEAT data is caused by one of those 
processes is not clear and can only be determined if fur- 
ther positron measurements provide improved statistical 
accuracy, and extend to higher energies. 

It has been suggested that at energies below a few GeV 
the measured positron fraction may reflect charge sign 
dependent effects of the solar modulation [iH and 
that discrepancies between results from different experi- 
ments can be related to this modulation effect. The most 
recent solar polar field reversal occurred just prior to the 
HEAT-pbar measurement and thus data from the three 
HEAT flights can be used to investigate this proposed 
effect. In the energy range 5-6 GeV - the lowest energy 
common to all observations - the measured positron frac- 
tion is consistent with what we measure in 1994 and 1995 
with HEAT-e± (see Fig. EJ. At 5 GeV, Clem et al. 
predict a decrease in the positron fraction by about 40%, 
which is inconsistent with our measurements. It would 
be of interest to reconcile our result with the measured 
dependence of the antiproton flux on solar activity by the 
BESS group |27|,,2a]. 

In conclusion, we find that a primary contribution to 
the positron intensity above a few GeV can still not be 
ruled out. 
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